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Capsule 

This paper describes a graduate level class project designed to develop understanding of aerosol 
optical properties and methods for their retrieval from satellite observations. 

Abstract 

Bridging the gap between current research and the classroom is a major challenge to today’s 
instructor, especially in the sciences where progress happens quickly. NASA Goddard Space Flight 
Center and the University of Maryland teamed up in designing a graduate class project intended to 
provide a hands-on introduction to the physical basis for the retrieval of aerosol properties from state-of- 
the-art MODIS observations. Students learned to recognize spectral signatures of atmospheric aerosols 
and to perform spectral inversions. They became acquainted with the operational MODIS aerosol 
retrieval algorithm over oceans, and methods for its evaluation, including comparisons with ground- 
based AERONET sun-photometer data. 
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1 . 


Introduction 


The direct and indirect radiative effects of aerosols on climate and climate change remain a 
source of uncertainty in climate research (IPCC, 2001). Aerosols (also known as particulate matter 
(PM)) play important roles in precipitation processes, reduced visibility, and human morbidity (Samet et 
al., 2000). Before the satellite era, information on aerosols came from limited surface based 
observations, which are not sufficient to describe their spatial and temporal variability of aerosols. 

With their vantage high above the earth, satellite observations of reflected and emitted radiances 
are increasingly being used to monitor aerosols and their interactions within the climate system (King 
and Kaufman, 2000). A highly successful project is the NASA’s Moderate Imaging Spectrometer 
(MODIS) (Salomonson et al., 1989), which has been observing from aboard the Terra (since 2000) and 
Aqua (since 2002) satellite platforms. By observing spectral radiances in 36 channels (from 0.412 to 
14.2 pm - http://modis.gsfc.nasa.gov) and at resolution ranging from 250 m to 1 km, MODIS is highly 
suitable for deriving aerosol properties over the oceans (Tanre et al., 1996; Tanre et al., 1997). The 
MODIS polar geo-synchronous orbit is such that most of the globe is covered each day, and completely 
covered every two days, thus providing comprehensive observations. Since the MODIS data are well 
characterized and easy to access, they are suitable for introducing students to the physics of remote 
sensing. 

As part of a graduate course in Remote Sensing, students used MODIS data and focused on the 
following: 

• Understanding how aerosols interact with reflected solar radiation. 
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• Obtaining working knowledge of inversion techniques for retrieving aerosol properties 
over oceans. 

• Implementing the MODIS algorithm to retrieve aerosol properties from observations. 

• Contrasting aerosol properties from different sites, and to relate them to geography and 

season. 

• . Evaluating satellite retrievals in the context of ground based aerosol measurements. 

This paper reviews the basic physics of remote sensing of aerosols and describes selected 
findings and lessons learned by the students. Data, codes, and detailed instructions needed to perform 
the exercises are available at: 

http://www.atmos.umd.edu/~levy/MODIS_Aerosol_Project 


2. Remote sensing of aerosols over the ocean from MODIS and AERONET 

Aerosols are suspended liquids or solids in the atmosphere, ranging in size from a few molecules 
to tens of micrometers. Aerosols having radii of between about 0.1 pm and 20 pm are of primary 
interest, where 1.0 pm denotes separation between “fine” and “coarse” Modes. The aerosol optical 
depth (AOD or x) defines the integral of the light extinction by aerosols within an atmospheric column. 
Spectral dependence of x is related to the aerosol size distribution (Eck et al, 1999), and the Fine Mode 
Weighting (FMW or r\) describes the contribution of the fine aerosols to the total x at a specific 
wavelength (X). By convention, x and rj are reported for X=0.55 pm. 

Sun photometers have been in operation for many years and when calibrated properly, provide 
accurate measurements of spectral x (function of X, i.e. xj. By accounting for known effects of 
molecular scattering (Rayleigh scattering) and gas absorption, sun-photometer measurements of spectral 
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direct beam extinction are directly related to x x as expressed by the Beer-Bouger extinction law. The 
current standard for ground-based sunphotometry is the AERosol Robotic NET work (AERONET) 
(Holben et al., 1998), providing continuous daytime measurements of t% with accuracy of ±0.01 at most 
wavelengths (Remer et al., 2005; Pinker et al., 2001; Eck et al., 2003) at over a hundred sites globally. 
Through measurements of spectral diffuse sky radiation, AERONET also provides estimates of aerosol 
size distribution (Dubovik and King, 2000). 

Whitby (1978) showed that a series of log-normal size distributions accurately describe aerosol 
sizes. In most situations, aerosols that interact with solar radiation can be modeled as a weighted 
combination of two log-normal modes, one for the fine and one for the coarse (Wang and Gordon, 1994) 
Theoretically, using radiative transfer codes, one can pre-compute a look-up. table (LUT) of simulated 
satellite-observed reflectance (normalized radiance), arising from a variety of aerosol conditions for a 
subset of realistic geometries (satellite and sun angles) and total aerosol loadings x (Tanre et al., 1997). 
For MODIS, there are four choices of fine modes and five choices of coarse modes, simulated for 2304 
solar/satellite angle combinations and seven values of x indexed by 0.55 p,m (Levy et al., 2003; Remer et 
al., 2005). Note that the actual spectral reflectance for each mode includes contributions from the 
surface and the molecular atmosphere (Rayleigh scattering). 

MODIS over-ocean aerosol properties (known as ‘Level 2’ products) are retrieved at 10 km x 10 
km resolution, meaning that the spectral reflectance values used for the retrieval must be representative 
of the target area. Remer et al., (2005) describes how the high resolution (250 and 500 m resolution) 
MODIS spectral reflectance observations (known as ‘Level 1’ data) are aggregated, corrected for gas 
absorption, and selectively averaged (including glint, cloud, and sediment masking), in order to 
represent the 10 km box. The aerosol retrieval assumes that this ‘Level 2’ spectral reflectance arises 
from the radiative effects of a combination of fine and coarse modes within the atmospheric column. 
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The total x is set to fit the observed reflectance at one wavelength, and r| is calculated to match the 
spectral dependence of x x . Incorrect values of ti lead to larger spectral ‘fitting errors’. A separate 
solution is found for each of the twenty combinations of fine and coarse modes, so that the ’best’ 
solution is the one combination of fine (‘f ) and coarse (‘c’) modes with the least fitting error. The 
spectral AOD (xj and rj are results of the inversion. 

Ichoku et al. (2002) developed an algorithm to operationally compare MODIS (Level 2) and 
AERONET products, by addressing the issue of spatial and temporal matching. Remer et al. (2002, 
2005) used this approach in ‘validation’ of MODIS total x retrievals over the ocean on a global scale. 
Remer et al., (2005) also discussed retrievals of x% and how they relate to retrievals of aerosol size 
parameters. In our class project, we used a set of nine MODIS/AERONET co-locations in order to 
understand how the MODIS inversion works, as well as to study the difference between MODIS and 
AERONET retrievals of x%. 

3. The class project 

The aerosol inversion, as summarized in a few paragraphs above, is actually quite complicated. 
It was assumed that students could best understand inversion and the resulting products by performing a 
simplified inversion “by hand.” Students used MODIS (Level 2) spectral reflectance as observed over 
the ocean near coastal and island AERONET sites (Figure 1) dicing 2001 and 2002. Coincident 
AERONET-derived x% data were used to evaluate the MODIS-derived products. In all, nine pairs of 
collocated MODIS/AERONET data were used for the project (Table 1). 

The skills required for this project are familiarity with FORTRAN compiling and use of “Y 
versus X” plotting software. A Fortran source code and data files with instructions for downloading, 
compiling and running the programs, are provided to the students and are now publicly available (at 
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http://www.atmos.umd.edu/~levy/MODIS_AerQsol_Project). The input data sets were formatted for 
easy import into a spreadsheet program (such as MS-Excel). 

4 . Results 

4.1 Aerosols and Spectral Radiances 

The students were given Terra-MODIS (Level 2) spectral reflectances at 7 channels for each of 
eight co-locations (first eight columns in Table 1). Figure 2 presents plots of reflectance versus 
wavelength (spectral reflectance) for the eight MODIS observations. All locations show decreasing 
reflectance with wavelength, primarily attributed to lower oceanic reflectance and molecular (Rayleigh) 
scattering at longer wavelengths. 

They found that different spectral dependence was related to location. Both Ascension Island 
and Capo_Verde have “flatter” spectral dependence than the other sites. Both sites are in the Atlantic 
Ocean and both (especially Capo Verde) are often in the path of the African dust aerosol (Carlson and 
Caverly, 1977; Grousse et al., 1988; Prosper©, 2001). Dust aerosols show lower spectral dependence 
than other types of aerosols (Eck et al., 1999) because of the larger number of coarse aerosols that 
scatter in the longer wavelengths. Venise and COVE are urbanized areas, and are characterized by 
urban/industrial aerosols (Remer et al., 1998; Dubovik et al., 2002). Due to the larger fraction of the 
fine mode in urban aerosols, these sites display larger spectral dependence (Eck et al., 1999). 

4.2 Hand-Calculated Aerosol Optical Depth 

This exercise was designed to introduce students to the MODIS aerosol retrieval algorithm, by 
having them manually perform a MODIS-Iike inversion. A simplified look-up table was provided for a 
pair of fine and coarse modes at selected geometries (from the www). Also provided was the MODIS 
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spectral reflectance for Gotland (ninth column in Table 1). The students performed the following tasks, 
(analogous to those within the operational algorithm). 

To start, they interpolated the entire LUT to the specific solar/ surface/ satellite geometry 
measured by MODIS (Figure 3). This specific-angle LUT was used to fit the observed spectral 
reflectance. This fitting was done by ‘halving’ iteration on r|, coincidently with interpolation on indexed 
total x within the LUT. For the first iteration, they attempted to match the satellite reflectance by 
assuming that the aerosol was of coarse mode only (i.e. r)=0.0). This resulted in total x at (0.55 pm) of 
about 0.19, with large differences (‘fitting error’) between calculated and observed spectral reflectance 
(Figure 4). Assuming only fine mode aerosol (rpl.O) led to lower fitting errors. Trying half and half 
(ri=0.5) showed significant improvement, while ‘halving’ again (rp0.75) resulted in an even better 
match with observations. Following Tanre et al. (1996), the iteration was considered complete after ten 
iterations and when the simulated reflectance was within 3% of the observations (Figure 5). In the case 
at hand, convergence was achieved after three to four iterations. The accepted solution was estimated to 
be r|=0.79 and x=0.22 (at 0.55 pm), with fitting error e of about 2.5%, but depended on how the students 
chose to handle significant digits in their calculations. 

4.3 Global aerosol retrievals and comparison with sun photometer data 

How do aerosol properties vary globally? The students were given instructions for compiling 
and running a stand-alone version of the operational MODIS algorithm (also available from the listed 
www). Using MODIS angle and reflectance data described in section 4.1 as input, the aerosol optical 
properties were retrieved for each case. Figure 6 compares the spectral dependence of x^, retrieved from 
both MODIS and AERONET at eight sites. 
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The spectral dependence of x was seen to resemble the observed spectral dependence of the 
reflectance for each location. This reflects the characteristics of the aerosols near each site. As 
expected, the spectral dependence of t x is smaller (flatter) for dusty or mid-ocean sites (Ascension 
Island, Capo Verde), than for more polluted (urban/industrial or smoke aerosol) near coastlines (COVE, 
Venice). Africa shows a combination of small (presumably smoke) and larger sea salt aerosols. The 
retrieval of r) mirrors the spectral dependence of x^, with the larger particles over Ascension Island and 
Capo Verde, and the smaller particles over COVE and Venise. These results are consistent with what is 
known about aerosol characteristics in these regions. 

Figure 6 shows good visual agreements for Venise #1, Villefranche and COVE. These three 
sites are from regions known to be dominated by fine mode, aerosols and the spherical aerosol 
assumption (Mie theory) is appropriate. On the other hand, for regions dominated by larger non- 
spherical dust particles, the differences between the two products are considerable (e.g., at Ascension 
Island and Dahkla). This implies that the aerosol models employed in the MODIS algorithm are not 
sufficient for describing the ambient aerosol properties at these sites (e.g. Levy et al., 2003). 
Interestingly, the differences between the MODIS and AERONET retrievals for regions with large 
aerosols tend to increase with increasing wavelength, while those for the regions with small aerosols 
tend to increase with decreasing wavelength. 

5. Summary and Conclusions 

The objective of this project was to expose graduate students to current research in aerosol 
remote sensing. The students learned about the spectral dependence of aerosol scattering and how it is 
exploited in the MODIS aerosol retrieval algorithm to estimate aerosol optical properties. 


9 



Due to the importance of aerosols on atmospheric radiative processes and their climatic 
implications, the project is relevant to interdisciplinary problems. Moreover, MODIS represents the 
state of the art science and technology for remote sensing of several earth-atmosphere processes. 
MODIS (as well as other satellite) data are free, easily available, and increasingly useful in many 
applications, so there is a great need to learn how to process such data and interpret them in new and 
exciting ways. Projects such as this one provide opportunity for students and young scientists to become 
familiar with (and less apprehensive of) datasets of this magnitude. 
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List of Tables 


Table 1 . MODIS/AERONET co-locations used for this project. The first eight (normal font) are 
used in the MODIS/AERONET comparisons. The ninth site (in italic) was used for the 
manual inversion. 
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List of Figures 

Figure 1 . Map of AERONET sites selected for this study . 

Figure 2. Plot of reflectance vs. wavelength (spectral reflectance) for eight MODIS observations. 
Data used for the later “manual” retrieval (Gotland) are plotted with squares. 

Figure 3. LUT table interpolated to MODIS observed Geometry. Filled circles are the modeled 

reflectance from the small mode for each of five indexed optical depths (0.2, 0.5, 1.0, 2.0 
and 3.0 at 0.55 pm), while the open squares are the modeled reflectance arising from the 
coarse mode. The filled triangle (dashed line) is the modeled Rayleigh reflectance, and 
the filled square (dotted line) are the MODIS observations for Gotland. Reflectance from 
each mode increases monotonically with x (AOD). 

Figure 4. Reflectance calculated “by hand” (colored circles) versus the reflectance observed by 
MODIS (black squares). Different colors represent different choices of the small mode 
weighting, r). The best fit to MODIS data (black squares) arises from r| = 0.787 (orange). 

Figure 5. Fitting error (green), fine mode weighting (red) and aerosol optical depth (blue) 
calculated at each iteration of the MODIS retrieval. 

Figure 6. Spectral ^retrieved by MODIS (solid) and AERONET (dashed) for each location. 
MODIS retrieved rj and fitting error are indicated on each plot. 
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Table 1. 


Date 

4/11/2001 

8/15/2002 

8/2/2001 

8/19/2002 

7/11/2002 

3/23/2001 

9/13/2002 

8/17/2001 

1/1/02 

DOY 

102 

227 

215 

231 

192 

83 

256 

230 

190 

Time (UTC) 

15:45 

11:55 

11:10 

11:15 

9:30 

11:45 

11:15 

15:25 

11:20 

Location 

Arica 

Ascension_ 

Island 

Venise#! 

Villefranche 

Venise #2 

CapoJVerde 

Dahkla 

COVE 

Gotland 


16 















M 

m _ 

m o 


Figure 1. 












Figure 3. LUT table interpolated to MODIS observed Geometry. Filled circles are the modeled ' 
reflectance from the small mode for each of five indexed optical depths (0.2, 0.5, 1 .0, 2.0 
and 3.0 at 0.55 pm), while the open squares are the modeled reflectance arising from the 
coarse mode. The filled triangle (dashed line) is the modeled Rayleigh reflectance, and 
the filled square (dotted line) are the MODIS observations for Gotland. Reflectance from 
each mode increases monotonically witht(AOD). 
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Reflectance calculated “by band” (colored circles) versus the reflectance observed by 
MODIS (black squares). Different colors represent different choices of the small mode 
weighting, % The best fit to MODIS data (black squares) arises from T] - 0.787 (orange) 
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